Equilibration and relaxation times at the chiral phase transition including reheating 
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We investigate the relaxational dynamics of the order parameter of chiral symmetry breaking, the 
sigma mean-field, with a heat bath consisting of quarks and antiquarks. A semiclassical stochastic 
Langevin equation of motion is obtained from the linear sigma model with constituent quarks. The 
equilibration of the system is studied for a first order phase transition and a critical point, where 
a different behavior is found. At the first order phase transition we observe the phase coexistence 
and at a critical point the phenomenon of critical slowing down with large relaxation times. We 
go beyond existing Langevin studies and include reheating of the heat bath by determining the 
energy dissipation during the relaxational process. The energy of the entire system is conserved. In 
a critical point scenario we again observe critical slowing down. 

PACS numbers: 



I. INTRODUCTION 

In heavy-ion collisions matter is created under extreme 
conditions, in small systems and with fast dynamics. It 
is, therefore, likely that nonequilibrium effects play an 
important role in the evolution of the fireball. Such ef- 
fects are of particular importance at the conjectured crit- 
ical point in the QCD phase diagram. Signals of the 
critical point are based on the growth of the correlation, 
which leads to large fluctuations [11, 01 ■ Naturally the size 
of the system itself is a limit to the correlation length. 
Here, the renormalization group methods of finite size 
scaling could be applied to localize the critical point Q. 
Finite time effects due to the phenomena of critical slow- 
ing down turn out to be more limiting on the growth of 
the correlation length than finite size effects [3j ■ The re- 
laxation times at a critical point become infinitely long. 
Even if the system is in thermal equilibrium above Tc it 
is necessarily driven out of equilibrium when crossing the 
critical point. 

For the dynamics in the vicinity of a phase transition 
or rapid crossover the study of an order parameter is of 
particular interest. For a proper description of nonequi- 
librium effects the following ingredients are necessary: 

• a (semiclassical) equation of motion for the long- 
wavelength modes related to the order parame- 
ter, which includes damping and noise for the 
long-wavelength modes, i.e. a Langevin dynamics 
needs to be formulated (see e.g. [1] and references 
therein).; 

• a fluid dynamic description of the rest of the system 
("heat bath") [1,0]; 

• the back reaction of the long-wavelength modes on 
the heat bath should be consistently included. 



This latter aspect is important because for a fireball with 
an extension of at most a few femtometers the "heat 
bath" cannot be regarded as infinitely large. But it 
is missing in most of the Langevin descriptions 
while chiral fluid dynamic models typically disregard the 
Langevin dynamics [1, [3] ■ 

In [l^ we have derived a Langevin-type equation of 
motion for a sigma field coupled to a heat bath of quarks. 
We have determined explicit expressions for the damping 
coefficient and the correlation of the noise term. Within 
the formalism of the two-particle irreducible (2PI) effec- 
tive action [l3l - [2l| we were able to obtain the equilib- 
rium properties of the heat bath consistently. During 
the relaxational process of the sigma field energy dissi- 
pates from the field to the heat bath, which is thereby 
reheated. Usually this reheating effect is not taken into 
account. In a small and dynamic system this back re- 
action on the heat bath can, however, be important and 
change the overall evolution of the system. 

In this paper we want to investigate the relaxational 
dynamics of the sigma field and time scales associated 
with it. The relaxational properties of the sigma field 
are important for studies of the phase transition. In a 
nonequilibrium situation one expects to see the coexis- 
tence phase for a first order phase transition, and due 
to critical slowing down one also expects long relaxation 
times at the critical point. 

In section |ll] we explain the dynamics for the sigma 
mean-field and the local equilibrium properties of the 
quarks on the basis of the linear sigma model with con- 
stituent quarks. We study the equilibration properties of 
the sigma field within our model for a static and isother- 
mal heat bath in section IIIII During the relaxation pro- 
cess the sigma field loses energy due to dissipation. This 
is studied in comparison with the energy conservation of 
the entire system in section IIVI Finally, we include this 
energy exchange and take reheating of the heat bath into 
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account in section [V] 

II. THE DYNAMICS OF THE SIGMA FIELD 

In our approach jl2| the relaxational dynamics of the 
order parameter of chiral symmetry breaking, the sigma 
mean-field, is given by the hnear sigma model with con- 
stituent quarks [l^. The quark degrees of freedom are 
assumed to be in local thermal equilibrium, such that 
they constitute the heat bath. The sigma mean-field is 
propagated by a stochastic Langevin equation of motion. 

A. The linear sigma model with constituent quarks 

The linear sigma model with constituent quarks de- 
scribes the chiral phase transition. Here, the a and tt 
mesons couple to quarks. The Lagrangian reads 

H = q[i^^d^, - g{a + 275x7?)] g 

+ \{d^ad^a) + \{d^T:d^T:) - t/(a,Tf) , 

with the constituent quark field q — {u, d), the sigma field 
a and the pion fields tt . g is the coupling strength of the 
interaction between the chiral fields and the quarks, f are 
the isospin Pauli matrices. The classical chiral potential 
is given by 

Uia,7:)^^{a^+7:'-,y'f-h,a-Uo. (2) 

The Lagrangian ([T]) is invariant under transformations 
of the chiral group SUl(2) x SUr(2) if the explicit sym- 
metry breaking term hq vanishes in the potential 
The parameters in ([2|) are chosen such that the vacuum 
expectation values are reproduced, where (cr) — fn — 
93 MeV and (7?) = 0. Here, chiral symmetry is sponta- 
neously broken. The explicit symmetry breaking term 
is hq — f-nTn^ with the pion mass 771-^ = 138 McV. 
Then, = — mJ/A^. A realistic vacuum sigma mass 
ml = 2\^fl + ml w 604 MeV is given for = 20. In 
order to have zero potential energy in the ground state 
the term Uq — to^/(4A^) — f'^ml is subtracted. At a 
coupling g = 3.3 the constituent quark mass in vacuum 
is mq = 306.9 MeV. 

The Lagrangian ([T]) treats the quarks and antiquarks 
and the mesons on equal footing. In the real world con- 
fining forces recombine quarks and antiquarks in mesons 
and baryons below the confinement critical temperature. 
The aspect of confinement is not included in the linear 
sigma model with constituent quarks. We can, thus, in- 
vestigate the pure effect of the chiral phase transition. 
For extensions of the modelgluons can be included on 
the level of the dilaton field U or the Polyakov loop [l^ . 

The linear sigma model with constituent quarks ex- 
hibits the full spectrum of the suggested chiral phase 
structure of QCD 24]. It has a crossover transition at 
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FIG. 1: The effective potential shown for the first order phase 
transition at a coupling g — 5.5 and the transition temper- 
ature Tc — 23.27 MeV and the critical point at a coupling 
g = 3.3 and the transition temperature Tc = 139.88 MeV. At 
the first order phase transition the effective potential has two 
minima, one close to vanishing a and one close to the vacuum 
expectation value (a) = fn- At the critical point the mini- 
mum is very flat around the minimum. The effective potential 
is given by the thermodynamic potential to be calculated in 

vanishing fi b and a first order transition line at high 
and lower t emp eratures, which terminates in a critical 
point at Tc ^M,- The phase transition can also be tuned 
by fixing /zs = and changing the coupling constant 
g [2^, [23] . At large couplings the phase transition is of 
first order. Here, the effective potential has two minima, 
which are separated by a barrier. Lowering the coupling 
constant one finds a critical point, where the curvature 
at the minimum of the effective potential is very fiat. In 
this work we use g — 5.5 for a first order phase transi- 
tion with the transition temperature Tc — 123.27 MeV 
and g = 3.63 for a critical point with the transition tem- 
perature Tc — 139.88 MeV. For both cases the effective 
potential is shown in figure [TJ The purpose of the present 
work is a first qualitative assessment of the relaxational 
dynamics of the order parameter consistently coupled to 
a heat bath of finite size. Therefore we did neither in- 
clude finite baryo-chemical potential nor pion degrees of 
freedom. Consequently we have to live with the some- 
what unrealistic values of the transition temperature. 

B. Local equilibrium properties of the quarks 

Within the formalism of the 2PI effective action p^| - [2]| 
the local thermodynamics of the quarks is given consis- 
tently by the one-loop thermodynamic potential in mean- 
field approximation. At /i^ = the thermodynamic po- 
tential is 

n{T, a,Tf)^U {a, tt) + nqq{T, a, 7?) 

= U (a, n) 2dqT /(^ In (1 + cxp 

(3) 
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where = 12 is the degeneracy factor of the quarks for 
Nf = 2 flavors, Nc = 3 colors and the two spin states. 
The quark mass is generated by nonvanishing expectation 
values of the chiral fields due to spontaneous symmetry 
breaking. In the course of the evaluation of the functional 
determinant in Dirac and isospin space one generates a 
term defined as the effective mass of the quarks 



m„ 



g [a +TT ) . 



(4) 



Then, the energy of the quarks and antiquarks is 

E = + m2 = y/]? + g^{a^ + tt^) . (5) 

According to the local value of the temperature and the 
sigma mean-field the local pressure of the quarks is 

p(cr,7r,T) = -Ogg(r,cr,7r) . (6) 

And the local energy density of the heat bath of quarks 



dp{a, 7?, T) 
e(cr,7r,r) =T — p[a,n,T) 



=24 



(2^ 



(7) 



Enp [p] 



C. Langevin equation of the sigma field 

The Langevin equation of motion of the sigma field 
can consistently be derived in the 2PI effective action 
formalism or, alternatively, in the infiuence functional 
method [2^. It gives a reduced description of the 
entire system with focus on the evolution of the relevant 
variables, which are propagated explicitly. The details 
of the environment are eliminated by integrating out the 
environmental fields in a path integral over the closed 
time path contour [29 . | 3C| . This method was widely used 
Isil - fssl. in gauge theories (s^, IsBj 
and in 0{N) chiral models [10] . 

For the linear sigma model with constituent quarks we 
assume the following splitting: the irrelevant degrees of 
freedom are the quarks and antiquarks, which constitute 
the heat bath, and the relevant sector is that of the sigma 
field, which we propagate explicitly. 

The Langevin equation for the sigma mean-field reads 

d^d^a{t, x) + ^ + ^ + r,{T)dta{t, x) = ^{t, x) . (8) 

OCT 0(7 

We calculated the damping coefficient and the noise cor- 
relator in 

For the kinematic range, m^lT) > 2mq(T) = 
2gacq{T), where the decay of the sigma in a quark- 
antiquark pair is allowed the damping coefficient rj for 
the zero mode of the sigma mean-field reads 



3/2 



(9) 



For the present calculations we use the respective 
equilibrium values of the sigma mass m^^^(T) = 

{d'^V,{T, a) / da^)\cr=a^q and crcq{T) to evaluate this crite- 
rion. We use ^ also as a approximation for the non-zero 
modes in the Langevin equation 

In this framework there would be no damping in the 
low-temperature phase because the quarks are not light 
enough to allow for the decay a qq. Physically this 
makes sense, because the quarks should be confined any- 
way. Additional damping is provided by the decay cr — > 
2TT. Strictly speaking this is not included in our present 
framework. Nonetheless, to obtain a more realistic setup 
we include the zero-temperature damping coefficient from 
[lo| when the constituent quark mass is too large for the 
a ^ qq decay. Then for ma{T) < 2mq{T) ~ 2gacq{T) 



rj ~ 3 /fm . 



(10) 



The stochastic field in the Langevin equation (jS]) has 
a vanishing expectation value 



(11) 



and the noise correlation is given by the dissipation- 
fluctuation theorem 

mm)i = ^^{t - t')m,VCoth (^) . (12) 

In this paper we investigate the time evolution of the 
following quantities. The volume average of the sigma 
field for one configuration of the noise field S^ijk is 



(o')ri — ^3 ^ 0-ijk,n , 



(13) 



ijk 



where N is the number of cells in each direction. We 
average over several different configurations of the noise 



(14) 



typically between Nr = 5 and Nr = 20 depending on 
how different the trajectories really are for the various 
temperatures. 

For the numerical implementation of the Langevin 
equation we apply the well tested algorithm used in [ll| . 
In the calculations presented here, we choose iV = 32. 
The size of time steps is At — 0.02 fm and the lattice 
spacing is Ax = 0.2 fm. 



III. EQUILIBRATION FOR A GLOBAL, 
ISOTHERMAL HEAT BATH 



First, we study the equilibration of the sigma field with 
a global, i.e space-homogeneous, and isothermal heat 
bath. Concerning the energy exchange the back reaction 
of the sigma field on the heat bath is ignored. There- 
fore, the temperature of the heat bath is constant and 
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determines the shape of the effective potentiaL This is 
very different for the first order phase transition and the 
critical point. We, therefore, expect a different evolution 
of the sigma field for these two scenarios. 

While the initial conditions for the sigma field are var- 
ied in the next sections, the time derivative of the sigma 
field dfd is initially zero. There is no clear physical mo- 
tivation for the choice of the initial dtcr. When it is ini- 
tialized in direction of the relaxation process the sigma 
field relaxes faster. A random distribution for the initial 
dt<J averages out. 




o-/MeV 



A. First order phase transition 

The evolution in a first order scenario is especially in- 
teresting because the effective potential has two minima 
in the spinodal region 108 MeV ~ T^l^ < < Tip^ ~ 
128 MeV, see figure [T] for the effective potential at the 
transition temperature Tc — 123.27 MeV. We expect that 
for some configurations the sigma field relaxes partly to 
the unstable minimum instead of the true thermal expec- 
tation value. This is even more likely in the vicinity of 
Tc, where both minima are almost degenerate. 

In order to study the relaxation of the sigma field to its 
thermal equilibrium state at the temperature of the heat 
bath we need a clear nonequilibrium initial situation. In 
our investigation we distinguish between two cases for 
which the initial nonequilibrium situation is realized dif- 
ferently: the equilibration at temperatures above and at 
temperatures below the transition temperature. 



FIG. 2: The distribution of the sigma field for T = 160 MeV 
at different times. The initial distribution is fiat. It quickly 
turns into a Gaussian centered at a mean of 2.89 MeV with a 
width of 1.62 MeV. The Gaussian fit is to the distribution at 
t = 80 fm 
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1. Equilibration at high temperatures T > Tc 

For the equilibration at temperatures above the phase 
transition temperature we initially distribute the sigma 
field linearly between (imin — and (Tmax > such that 
the initial average of the sigma field is ((t)„ ~ 50 MeV. 
The flat distribution is far from being thermal. 

In figure [2] we see how the initially fiat distribution 
develops for an example evolution at T = 160 MeV. It 
becomes Gaussian after times t ~ 20 fm, which corre- 
spond to the relaxation times in figure |31 It shows the 
equilibration of the sigma field to its proper equilibrium 
value at the respective temperature. This is ensured by 
the fluctuation-dissipation theorem ((T^), which relates 
the variance of the noise field to the damping coefficient. 

The dynamics of the system becomes different for 
Tsp > T > Tc- In figure[4]the scaled noise average (ct) /coq 
and the volume average (cr)n/(Tcq for two individual noise 
realizations are shown for T — 125 MeV. Because of the 
two minima the system takes a long time to relax to the 
global minimum at tXeq — 6 MeV. Since the initial average 
is slightly shifted towards the low-temperature minimum 
the system first tends to this phase, but finally relaxes to 
the true minimum at very large times. 



FIG. 3: The time evolution of the scaled noise average of the 
sigma field in a first order scenario for temperatures above 
the upper spinodal temperature T > Tap . 
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FIG. 4: The long relaxation proces s of the sigma field for T — 
125 MeV. The scaled noise average (o")/o"cq and two individual 
noise configurations are given. Due to the barrier between the 
two minima the relaxation to the global minimum is very slow. 
The individual noise configurations differ substantially. 
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FIG. 5: The distribution of the sigma field at the transition 
temperature Tc = 123.27 MeV of the first order phase transi- 
tion scenario. The system shows phase coexistence. 
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FIG. 6: The time evolution of the scaled noise average of the 
sigma field in different quenched scenarios for a first order 
phase transition scenario. 



At the transition temperature Tc = 123.27 MeV the 
distribution of the sigma field for one noise configura- 
tion is shown in figure [5] The system is in the expected 
phase coexistence and the sigma field does not relax at 
all. Instead we observe that the system is split into one 
part in the high-temperature and one part in the low- 
temperature minimum. 



2. Equilibration for a quench to temperatures T < Tc 

After initializing the sigma field in equilibrium with 
an initial temperature of Tini = 160 MeV > Tc — 
123.27 MeV the system is suddenly quenched to different 
temperatures T < T^ for which the Langevin equation 
for the sigma field is solved. 

It is known that for the linear sigma model with con- 
stituent quarks the nucleation rates are rather low and 
that the main relaxation mechanism is that of spinodal 
decomposition ^26 ] . This is also observed in our calcula- 
tions. In figure [6] we show the time evolution of the re- 
laxation of the sigma field. For T = 115 MeV, where we 
still have a substantial barrier between the two minima, 
the relaxation times are significantly larger than for lower 
temperatures. Here, the system remains in the local min- 
imum cr ~ 10 MeV until at i ~ 15 fm it begins to decay 
to the global minimum, which is also a slow process. It 
resembles the case of an exponentially damped system, 
which decays without oscillations. At T = 100 MeV we 
can clearly see the oscillating relaxation. It occurs much 
faster. 



B. Critical point 

The effective potential for a scenario with a critical 
point has only one minimum at all temperatures. It 
continuously shifts from cr ~ to the vacuum expec- 
tation value (tT)vac = /tt when lowering the temperature. 
At the critical point this minimum becomes very flat. 
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FIG. 7: The time evolution of the scaled noise average of 
the sigma field for a critical point scenario for various tem- 
peratures. The relax;ation time becomes large at the critical 
point. 

see figure [1] and we expect long relaxation times. For 
different temperatures, the evolution of the noise aver- 
aged sigma field can be best compared by choosing the 
same initial conditions, the thermal equilibrium state at 
T = 160 MeV. The system is then quenched to temper- 
atures T < 160 MeV. The results are shown in figure 
[T] For low temperatures the potential is steeper and the 
relaxation process occurs faster. The field oscillates be- 
fore relaxing. Approaching the critical temperature re- 
laxation times become larger with a clear maximum at 



IV. THE ENERGY DISSIPATION 

During the relaxation of the sigma mean-field to its 
equilibrium value the dissipative term in the Langevin 
equation ([5]) causes energy dissipation. By the interac- 
tion with the quarks it is transfered to the heat bath. In 
[l^ we derived a conserved energy-momentum tensor of 
the entire system including a dynamics of the heat bath. 
In this paper, we do not include the fiuid dynamic ex- 



6 



pansion of the heat bath. In this section we investigate 
the relevant energy exchange between the sigma field and 
the heat bath. The energy dissipation of the field to the 
heat bath can be obtained from the energy-momentum 
tensor of the sigma field 



T]dt(j)dt(J, 



(15) 



where T^'' is the energy-momentum tensor of the purely 
mesonic Lagrangian 



Then, the energy dissipation is described by 

AE^diss = {gps + r]dt<j)dt<jAt . 
The total energy of the sigma field is given by 



1 



U{a) 



(16) 



(17) 



(18) 



It has a kinetic, potential and fluctuation energy term. 
During relaxation to the vacuum expectation value the 
potential energy is transfered to kinetic energy as dtcr 
grows. Then, the damping becomes substantial and 
causes energy dissipation. This flow of energy from the 
fleld to the heat bath is given by ([TT]) . There is a reverse 
flow of energy from the heat bath to the field AE^ as- 
sociated with the noise field ^ in the Langevin equation 
(|S]), which is an averaged quantity balancing the energy 
dissipation AE'diss in equilibrium and thus restoring the 
proper thermal equilibrium. This was already discussed 
in (l2j . Assuming that the made approximations in [l2| 
cause only a small violation of energy conservation we can 
determine A£'^ from comparing AE'diss to the energy dif- 
ference in the fleld before and after each numerical time 
step AEa- 

Here, we first show that the difference between Ai?diss 
and AE^ is small if one ignores the noise term in the 
Langevin equation. This is shown in figure 8(a) for the 
quench from Tjni = 160 MeV to T = 100 MeV in a sce- 
nario with a first order phase transition and for a crit- 
ical point scenario quenched from Tini = 160 MeV to 
T = 130 MeV in figure [8(b)] The resulting difference 
is a measure of the violation of energy conservation due 
to the approximations made in [l2j. It is numerically 
small. In figure |5] we also show the difference between 
Ai?diss and AE^^ including the noise term. We identify 
this difference with AE^. 



EQUILIBRATION FOR A HEAT BATH 
WITH REHEATING 



In this section we investigate the influence of the en- 
ergy conservation on the equilibration of the entire sys- 
tem. While the sigma fleld relaxes after a sudden tem- 
perature quench energy dissipates from the system to the 
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(b) 



FIG. 8: Energy diss ipation for a scenario with a f irst order 
phase transition (a) and for a critical point scenario (b) The 
system is quenched from Tini = 160 MeV to T = 100 MeV 
in the first order phase transition scenario and from Tini = 
160 MeV to T = 130 MeV in the critical point scenario. The 
Langevin equation is once solved with the noise term ^ and 
once without it. For each case the comparison between the 
total energy dissipation "Y^y Ai?diss and the energy difference 
in the field AE^ summed over the whole volume is shown. 



heat bath. This in return changes the temperature of the 
quark fluid and the effective potential. Thus, the evolu- 
tion of the sigma fleld itself is altered. In the last section 
we discussed the energy transfer between the sigma field 
and the heat bath. It has the two components Ai^diss 
and AE^. In the following we locally calculate AE^r and 
add this to the local energy density of the heat bath given 
by ([7]). The new energy density is inverted to flnd the 
local temperature. 



A. First order phase transition 

We present four results for scenarios with a flrst order 
phase transition. We quench from Tini — 160 MeV to 
Tsys = 100, 80, 50 and 20 MeV. During the relaxa tion 



of the volume averaged sigma fleld, see flgure |9(a)[ the 
average temperature increases rapidly to Tfin, see flgure 
10(a) The exact values are shown in table (a) Three 
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FIG. 9: Time evolution of the volume averaged sigma field 
in a scenario with a first order phase transition (a) and with 



a critial point (b) for different temperature quenches. The 
energy dissipation from the sigma field to the heat bath is 
taken into account. 



FIG. 10: Time evolution of the temperature in a scenario 



with a first order phase transition (a) and with a critical point 



(b) for different temperature quenches. The temperature is 
changed by the energy dissipation from the sigma field to the 
heat bath. 



temperatures are above Tc and below or close to the up- 
per spinodal temperature Tsp , where the effective po- 
tential has two minima. The sigma field initially relaxes 
towards the vacuum value. This relaxation reheats the 
heat bath and causes an increase in the temperature to 
above Tc. Large parts of the sigma field now remain in 
the unstable low-temperature minimum. We see that in- 
cluding reheating the entire system does not equilibrate 
for these temperatures. Obviously, the reheating locally 
changes the effective potential such that it counteracts 
the relaxational process. Only for the very low tempera- 
ture Tsys — 20 MeV, which is close to vacuum conditions, 
the final temperature is below T,.. Thus, the initial re- 
laxation of the sigma field corresponds already to the 
equilibrium state at Tfin. 



B. Critical point 

In a scenario with a critical point the effective poten- 
tial has only one minimum for all temperatures. There- 
fore, we expect the entire system to equilibrate. We 
consider the following four temperature quenches from 



160 MeV to Ts, 



140, 130, 100, and 80 MeV re- 



spectively. The corresponding volume averaged values 
are shown in table |(b)[ The volume averaged variances 
of both quantities are explicitly given. We clearly see 
that the entire system relaxes at a temperature Tf^^^ and 
(T = Tfin) in figures [9(b)] and [TO(b)| We ob- 
serve that for a temperature quench to Tgys = 130 MeV 
the final temperature comes closest to the critical tem- 
perature Tc = 139.88 MeV. As seen in figure 9(b) and 
figure [lO(b)[ relaxation times are longest for this quench. 



VI. SUMMARY AND OUTLOOK 

We have studied the dynamics of the sigma field given 
by the Langevin equation ([5]) with the damping ^ and 
(|T0)) and the noise correlator (IT2]) . It leads to the re- 
laxation of the sigma field with a static isothermal heat 
bath. Including reheating of the heat bath we find full 
relaxational dynamics only for a scenario with a critical 
point. For a first order phase transition the system stays 
in the low-temperature minimum, which due to reheating 
becomes the unstable minimum at the final temperature. 
At the critical point we observed longest relaxation times 
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Tfln/MeV 


rjfin/MeV 


<Tcq(r = rfln)/MeV 


100 


131.89 


69.06 


5.00 


80 


127.69 


80.11 


5.56 


50 


123.59 


87.57 


6.25 


20 


121.41 


85.60 


84.82 



(b) 



T,y,/MeV 


Tfln/MeV 


(Tfln/MeV 


aeq(T = Tfl,)/MeV 


140 


141.42 ±0.31 


23.33 ± 1.92 


24.37 


130 


138.96 ± 0.47 


53.43 ± 1.52 


54.52 


100 


124.53 ±0.57 


78.46 ± 1.40 


78.60 


80 


115.44 ±0.62 


83.82 ± 1.43 


83.90 



which is in accordance with critical slowing down. In fu- 
ture work we will include the fluid dynamic expansion 
of the heat bath and thus study the full nonequilibrium 
chiral fluid dynamics. 



TABLE I: Exact values for the relaxation of the volume aver- 
aged sigma field and the final temperatures for the different 
quenches for Tini = 160 MeV to Tsys. Here for a scenario with 



a first order phase transition in (a) and with a critical point 
in [(b)] 
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